Methyl group dynamics in a confined glass 
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We present a neutron scattering investigation on methyl group dynamics in glassy toluene confined 
in mesoporous silicates of difi'erent pore sizes. The experimental results have been analysed in 
terms of a barrier distribution model, such a distribution following from the structural disorder in 
the glassy state. Confinement results in a strong decreasing of the average rotational barrier in 
comparison to the bulk state. We have roughly separated the distribution for the confined state in 
a bulk-like and a surface-like contribution, corresponding to rotors at a distance from the pore wall 
respectively larger and smaller than the spatial range of the interactions which contribute to the 
rotational potential for the methyl groups. We have estimated a distance of 7 A as a lower limit of 
the interaction range, beyond the typical nearest-neighbour distance between centers-of-mass (4.7 

A). 

61.43.-j Disordered solids, 61.12.-q Neutron diffraction and scattering, 61.46.-fw Nanoscale materi- 
als: clusters, nanoparticles, nanotubes and nanocrystals 



I. INTRODUCTION 

Porous materials have attracted great interest in the 
last years, since they provide a framework for the study 
of the physical propierties of solids under restricted 
geometry^. The dynamics of systems confined in these 
materials is expected to differ considerably from that in 
the bulk state, first of all because the length of the inter- 
molecular interactions and correlations is limited by the 
finite size of the pore if the latter is sufficiently small. 
Moreover, the interactions between the pore wall and the 
confined molecules -which strongly depend on the porous 
system- would affect the dynamics of these molecules in 
a non-well defined surface layer on the confined system. 
Due to their simple dynamics^, small molecular rotors as 
ammonia, methane, ammine ions or methyl groups can be 
considered as adequate dynamical probes of the effects of 
confinement by comparing the rotational barriers in the 
bulk and in confined geometry'^. 

A thorough neutron scattering study on methyl groups 
dynamics in bulk glassy toluene has been recently pre- 
sented by some of US'*. The experimental results were 
successfully analysed in terms of a well-established bar- 
rier distribution model, such a distribution resulting from 
the structural disorder present in the glass. The aver- 
age barrier of the distribution was found to be notably 
higher than that of the /3-crystalline reference phase, 
which shows a short-range structure similar to that of 
the glassy phase"'"^. This fact could indicate that the 
relevant intermolecular interactions contributing to the 
rotational potential of the methyl groups extend beyond 



nearest- neighbours . 

In order to make a comparison with the bulk state and 
to shed new light about this question and about the in- 
teraction lengths controlling the parameters of the distri- 
bution, we have carried out neutron scattering measure- 
ments on glassy toluene confined in mesoporous silicates 
with different pore sizes. 



II. THEORETICAL ASPECTS 

The usual single-particle model for methyl group dy- 
namics in crystalline systems at very low temperatures 
(T « 2 K) is that of a rigid rotor tunnelling through a 
one-dimensional rotational potential V^(*&), which is re- 
stricted to keep the rotational symmetry of the methyl 
group^. In most of cases only the leading threefold term 
of the Fourier expansion is retained, higher order cor- 
rections being small: V{'^) = ^3(1 - cos3$)/2. The 
torsional levels are tunnel-split by the overlapping of the 
neighbouring single- well wavefunctions, and the energy 
splitting of the ground torsional level, or rotational fre- 
quency Wt , can be observed in neutron scattering spectra 
-for barrier heights V3 below w 700 K- in the /ieV-range 
as two resolution- width inelastic peaks centered at itJt. 
The corresponding incoherent scattering function for ro- 
tational tunnelling, normalized to scattering from one hy- 
drogen is^: 



S'^'o{Q,uj) = [A{Q) + B{Qmoj) 
+B{Q)[5{u: + ujt) + 5{u: - ujt)] 



(1) 
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with hQ^ huj respectively the momentum and energy 
transfer of the neutron, A{Q) = [1 + 2jo(<9r)]/3, B{Q) = 
2[1 - jo{Qr)]/9, joiQr) = sm{Qr)/{Qr) and r the H-H 
distance in the methyl group. 

At temperatures typically above ~ 50 K, a picture 
of thermally activated classical hopping over the bar- 
rier is valid. The corresponding hopping rate fol- 
lows an Arrhenius-like temperature dependence F = 
T^exp{—EA/kT), with Ea the classical activation en- 
ergy, defined as the difference between the top of the bar- 
rier and the ground torsional state, r^o is a temperature- 
independent preexponential factor. The incoherent scat- 
tering function for the classical hopping regime is^: 



S'uGiQ,^) = A{Q)6{uj)+3B{Q)L{uj;T) 



(2) 



with L{uj:T) a normalized Lorentzian of half- width at 
half-maximum (HWHM) equal to F. 

A scries of investigations in structural glasses has ev- 
idenced the need of introducing a barrier distribution 
giVs) to give account for the different features observed in 
the spectra of these systems (rotation-rate-distribution- 
model, RRDM)^. Such a distribution has its origin in the 
different local environments for the methyl groups result- 
ing from the structural disorder in the glassy state. Ac- 
cording to this approximation, the spectra for the glassy 
system are obtained as a superposition of crystal-like 
spectra (tunnel-like as (1) at T w 2 K, or classical-like as 
(2) at high temperature) weighted by (/(VJj): 



S'^h{Q,UJ,V3)dV3 



(3) 



Numerical relationships^ between V3, hivt and Ea al- 
low the transformation between the corresponding dis- 
tributions of these quantities ^(Vs), h{fiLOt) and f(EA)- 
By assuming that the preexponential factor Too is V3- 
independent^'®, /(-Ea) can also be transformed, through 
the Arrhenius law for the temperature dependence of the 
hopping rate, into the corresponding distribution of hop- 
ping rates iJ(logF). Thus, the general equation (3) for 
the RRDM is reduced to the tunnelling and hopping cases 
by doing the substitutions 5{u)±uit) J h{±Jiuit)d{hu>t) 
in equation (1) and L{uj,T) -> J H(logF)i(u;,F)(i(logF) 
in equation (2). The consistency of the model is reached 
when the distributions of tunnelling frequencies and clas- 
sical hopping rates, obtained respectively at low and high 
temperature, follow from the same barrier distribution 



III. EXPERIMENTAL DETAILS 

Neutron scattering measurements were carried out at 
the backscattering spectrometer IN16 of the Institute 
Laue Langevin (ILL, Grenoble, France). A wavelength 
of 6.27 A was selected providing an energy resolution 
of 0.5 yueV (HWHM). The energy window available by 



Doppler shift covered from -15 to 15 /ieV. The scatter- 
ing angle covered a range from 11 to 149°, resulting in a 
Q-window from 0.2 to 1.9 A~^. The instrumental resolu- 
tion was calibrated by a vanadium sample, which shows 
purely elastic scattering. Raw data were corrected for de- 
tector efficiency, scattering from the sample holder and 
absorption by means of ILL standard programs. 

Spectra were taken on three samples of ring-deuterated 
toluene confined in mesoporous silicates of the series 
MCM-41 and SBA-15, which geometry consists of a 
honeycomb-type lattice of parallel cylindrical pores. 
They were synthesized with pore diameters D = 24, 35 
and 47 A, according to an hydrothermal procedure de- 
scribed in Ref. [7]. The porous geometry parameters 
were confirmed by neutron diffraction and nitrogen ab- 
sorption experiments. Deuteration of the silanol groups 
was achieved by chemical H/D exchange with no further 
chemical treatment of the surface. Complete filling of 
the outgased matrices was achieved with an appropriate 
mass of toluene^. 

Flat samples of thickness 2 mm were used to get trans- 
missions above 90 %, allowing to neglect multiple scat- 
tering effects in the Q-range investigated, 1.3 - 1.9 A , 
where the ratio B{Q)/A{Q) is most favourable in the 
Q-window of IN16. The cooling rate from room temper- 
ature to 2 K achieved in the cryostat was sufficient to 
get the glassy state of confined toluene. Spectra were 
also taken for the empty matrixes, showing purely elastic 
scattering in the investigated temperature interval. 



IV. RESULTS 

Figs. 1,2 show the incoherent scattering functions, nor- 
malized to maximum unity, for confined and bulk glassy 
toluene -the latter taken from the previuos investigation 
of Ref. [4]-, obtained after substraction of all the elastic 
contributions due to coherent scattering and to incoher- 
ent scattering from other atoms different from the methyl 
group protons. In this way the comparison between the 
dynamics in the confined and in the bulk state is not af- 
fected by the different ratios from the methyl group inco- 
herent cross section to the total one. From a direct visual 
inspection, it is evidenced that confinement results in a 
strong change in the tunnelling dynamics at 2 K (Fig. 1), 
and also at the hopping dynamics at high temperature 
(Fig. 2). As previously reported for the bulk state*, the 
experimental data for the confined samples at 2 K could 
be reproduced in terms of a superposition of tunnelling 
frequencies, as exposed in Section 2, such a distribution 
following from an asymmetric Gamma-distribution of ro- 
tational barriers ^(Vs): 
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with r(p) the Euler Gamma function, ¥3^ the average 
barrier and p an adimensional shape parameter measur- 
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ing the asymmetry of the fimction. The corresponding 
standard deviation cry is obtained as ay = V3„/^. The 
thick lines in Fig.l are tlie RRDM theoretical functions 
for rotational tunnelling after convolution with the in- 
strumental resolution. The obtained parameters for the 
three pore diameters investigated are given in Table 1. 
The error bars are ~ 5%. The resulting distributions are 
of similar width, but much more asymmetric and with 
a notably lower average barrier than that of the bulk. 
A systematic increase of the average barrier is obtained 
with growing pore size, though far from the bulk limit in 
all cases, suggesting an important contribution of surface 
effects to the barrier distribution. 

In order to check the consistency of the RRDM for the 
obtained parameters of giVs), we derived the distribution 
of classical hopping rates through Arrhenius relationship 
from the distribution of classical activation energies de- 
rived from 5(^3). The corresponding RRDM theoretical 
function for classical hopping was constructed for differ- 
ent values of Fqc , convoluted with the instrumental reso- 
lution and compared with the experimental data. Thick 
lines in Fig. 2 correspond to the values Foo = 9 and 
8 meV for confined and for bulk glassy toluene^ respec- 
tively. The error bar for can be estimated in the 
order of ± 3 and ± 2 meV for the confined and the bulk 
cases respectively. 

V. DISCUSSION 

The single-particle potential for methyl group rota- 
tion in a given system is built up by intra- and inter- 
molecular interactions extending until a certain radius 

Tc. The results above described show that the distri- 
bution of potential barriers for methyl group rotation 
in glassy toluene is strongly modified by confinement. 
This has to be certainly related with the presence of 
the pore-wall and the interactions between the wall and 
the methyl groups. In a rough approximation we can 
expect that the methyl groups which are at a distance 
from the pore wall r < Tc will feel the wall and their 
single-particle potential will become modified. On the 
other hand, methyl groups at r > Tc from the pore wall 
will in principle behave as bulk-like methyl groups. In 
this framework, g{V^) can be separated into a "bulk-like" 
and a "surface-like" contribution, corresponding, respec- 
tively, to molecules at a distance of the pore surface be- 
yond and below Vc. Previous results corresponding to 
dynamics of small rotors in confined crystalline systems 
seem to support this kind of separation^. Prom the cylin- 
drical geometry of the pores, and assuming an homo- 
genous density of methyl groups within the pores, it is 
straightforward to calculate the fraction of surface-like 
molecules as /s = 1 — (1 — Ir^jDY ■ The corresponding 
surface-like contribution to the total barrier distribution 
can be calculated as.gs(V:3) = [,9(^3) - (1-./s).9b(V"3)]//s, 
with g{yz) and g^iy-i) the distributions for confined and 



bulk toluene defined by the parameters given in Table 1. 
For Tc below a certain value, this procedure yields to un- 
physical negative values of gs(V3), so that value can be 
understood as a lower limit of r^- Unphysical distribu- 
tions were obtained for Tc < 5, 6 and 7 A for D = 24, 35 
and 47 A respectively. Thereby these values correspond 
to the maximum possible bulk- like contribution to giy-i) 
in each case. Fig. 3 shows the distributions for one of 
the pore diameters and for the bulk case, and also the 
surface contributions for different values of r^, including 
some of those with unphysical negative values due to a 
underestimated value of (sec inset). As the interac- 
tion radius Tc should not depend on the pore size and 
Tc =7 A is also compatible with the g{y-i) obtained for 
the pore diameters = 24 and 35 A (see Table II), we 
can take 7 A as a first estimation for in glassy toluene. 
It is noteworhthy that this value has to be considered 
only as a lower limit, because larger values could also be 
compatible with the different giVz) obtained. Measure- 
ments at larger pore sizes would be required to give a 
more accurate estimation of rc. In any case our results 
seem to support the idea that the relevant interactions 
on the methyl group extend beyond the typical nearest- 
neighbour distance between centers-of-mass (4.7 A)^'^, 
as was suggested by the comparison of the results in the 
bulk glassy and /3-crystalline phases^. 
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TABLE I. 


Parameters 


of the barrier distribution for the 


different pore 


diameters and for the bulk case. 


D{k) 


Vf„(K) 


A (K) 


P fs 


24 


224 


155 


2.1 0.83 


35 


187 


124 


2.3 0.64 


47 


183 


119 


2.4 0.51 



TABLE IL Parameters of the surfacc-Iiko contribution 
5s(V3) calculated by assuming rc = 7 A. The fraction of sur- 
face-like methyl groups for this value rc is also included. 
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FIG. 2. As Fig. 1 for high temperature. In this case, thick 

lines are the theoretical functions given by the classical version 

of the RRDM. 



-0.015 -0.01 -0.005 0.005 0.01 0.015 
energy transfer (meV) 
FIG. 1. Incoherent scattering function for methyl group 
dynamics at 2 K for confined and glEissy bulk toluene. Cir- 
cles are the experimental data. Thick lines are the theoretical 
functions given by the tunnelling version of the RRDM convo- 
luted with the instrumental resolution (see text). Thin lines 
are the instrumental resolution. Amplification scale: 3% of 
the maximum. 
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FIG. 3. Barrier distribution .g(V^j) for confined glassy 
toluene with D = 47 A and its separation into a bulk-like 
and a surface- like contribution for Tc = 7 A (see text). The 
average barrier is indicated in all cases. The inset shows the 
surface-like distribution gs{Va) for different values of r^. The 
horizontal dashed line marks the zero value for gsiVs). 
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